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摘要 :以 天 然 草场 中 毒 杂 草 与 可 食 牧 草 为 研究 对 象 , 在 种 间 竞 争 的 生态 数学 模型 的 基础 上 引入 和 人 侵 扩散 因素 建立 毒 杂 草 入侵 扩 
散 模型 。 采 用 元 胞 自动 机 理论 将 竞争 模型 扩展 到 空间 网 络 进行 模拟 研究 ,分 析 毒 杂 草 属 的 空间 分 布 类 型 为 毒 杂 草 的 控制 提供 
数据 支持 。 研 究 表明 :(1) 在 入侵 扩散 作用 下 , 毒 杂 草 与 可 食 牧 革 的 共存 平衡 点 由 兽 加 为 两 个 味 兽 加 了 共存 的 可 能 性 ;(2) 
入 侵 扩散 作用 影响 了 毒 杂 草 种 群 的 空间 分 布 特征 ,减少 了 种 群 空间 分 布 的 聚集 程度 。 
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Abstract: Im the present study, a spatial diffusiom model of poisonous weed invasion by using inter-species competition 
model was developed based on poisonous weeds and edible forage in native pastures. Using the cellular automata theory, the 
competition model was extended to the sfatial\network to simulate and study the spatial distribution patterns of poisonous 
weeds, which will provide data supportihg their control. The results showed that (1) under the influence of invasion and 
diffusion, the coexistence eqtiilibftium point of the poisonous weeds and the edible grass could be increased by an increase of 
two, which increased the possibility of coexistence; and (2) the spatial distribution of the weed population was affected by 


invasion and diffusion, and the aggregation of the spatial distribution of the population was reduced. 
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生物 入 侵 是 指 生 物 由 原 分 布 区 域 经 由 自然 的 或 人 为 的 途径 扩展 到 一 个 新 的 区 域 ,繁殖 并 维持 下 去 '"。 
近年 来 , 随 着 全 球 经 济 一 体 化 .国内 外 贸易 的 迅速 发 展 , 世界 各 地 生物 入 侵 现象 愈加 严重 。 我 国 因 外 来 物 
种 入侵 造成 的 损失 相当 惊人 ,每 年 仅 几 种 外 来 物种 造成 的 经 济 损失 就 达到 574 亿 元 人 民 币 ”。 据 世界 自然 保 
护 联盟 报道 ,全 球 范 围 内 入 侵 物 种 对 生物 多 样 性 的 威胁 高 居 前 5 位 , 外 来 生物 入 侵 是 导致 原生 物种 衰竭 、 生 
物 多 样 性 减少 的 重要 原因 之 一 和 。 如 今 , 生物 入 侵 对 我 国生 物 多 样 性 的 影响 越 来 越 突出 ,在 一 定 程度 上 威 
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胁 我 国生 物 多 样 性 和 生态 系统 稳定 性 。 

我 国 西北 地 区 因 其 独特 的 地 形 条 件 和 地 理 位 置 , 加 之 资源 过 度 开发 和 不 合理 利用 ,导致 长 期 以 来 生物 灾 
害 严重 .生态 环境 脆弱 ,而 外 来 物种 入 侵 加 剧 了 这 一 了 矛盾。 外 来 和 人 侵 植物 可 降低 土壤 的 营养 水 平 , 这 主要 是 
由 于 竞争 .落叶 的 营养 贫乏 或 难 分 解 , 积累 盐分 改变 土壤 pH 值 等 造成 的 汪 。 在 天 然 草场 中 ,入 侵 的 毒 杂 草 
不 仅 与 可 食 牧 草 争夺 水 .肥料 .光照 等 资源 抑制 可 食 牧 草 的 生长 .降低 可 食 牧 草 的 产量 ,而 且 影响 可 食 牧 草 的 
品质 ,更 有 甚 者 ,有 的 杂 草 本 身 就 是 毒 草 » 家 冀 误 食 ) 造成 中 毒 .死亡 。 原 产 于 中 美洲 的 紫 葵 泽 兰 ( Eupatorium 
odoratum) , 约 在 20 世纪 50 年 代 初 从 中 缅 、 中 越 边境 传 入 云南 南部 ,至 今 已 占据 了 四 川 深山 州 175 的 草场 , 正 
以 每 年 约 30 km 的 速度 向 北 ,向 东 扩 散 草 延 ,当地 畜牧 业 遭 受 了 空前 打击 ,每 年 经 济 损失 达 数 千石 元 筷 。 丙 
此 ,牧民 都 希望 草场 中 的 毒 杂 草 越 少 越 好 ,可 食 牧草 越 多 越 好 ,而 要 想 控制 草场 中 毒 杂 草 的 数量 , 束 首 移 需 要 
了 解毒 杂 草 与 可 食 牧 草 在 草场 中 的 空间 分 布 情况 。 

本 文 假设 毒 杂 草 与 可 食 牧 革 存 在 对 环境 资源 的 竞争 关系 ,采用 元 胞 自动 机 理论 对 种 间 兑 争 模 型 进行 空间 
模拟 。 元 胞 自动 机 的 基本 思想 就 是 将 系统 分 解 为 多 个 元 胞 ,每 个 元 胞 看 成 一 个 单元 ,| 系统 中 各 个 相 邻 单元 之 
间 存在 一 定 的 交互 相互 关系 ,而 每 个 单元 的 状态 都 是 整个 系统 在 一 些 简单 的 规则 着 用 下 不 断 进 化 的 
结果 '*”。 

在 对 种 间 苋 争 模型 进行 模拟 的 时 修 , 本 文 引 入 了 空间 扩散 因素 ,从 而 得 到 了 毒 杂 和 草 和 人 侵 的 空间 分 布 情况 。 


1 毒 杂 草 与 可 食 牧草 的 种 间 竞 争 模型 


假设 毒 杂 草 与 可 食 牧草 种 间 存 在 相互 作用 ,同时 假设 单一 物种 独自 在 同一 个 草场 中 生存 时 ,数量 的 演变 
均 遵 从 Logistic 规律 。 于 是 得 到 毒 杂 草 与 可 食 牧 草 的 种 间 竞 争 模型 为 “” 


| ~ a 二 ) 

dt Ni N, 

和 (1) 
由 % ¥ 

一 《二 kk 

df ry NWN, Fl 


式 中 ,r， > 0,m > 0,a > 0,a，> 0, 记 Ww 分别 是 可 食 牧 革 和 毒 杂 草 的 数量 , rn , 7, 分 别 是 它们 的 固有 增长 
率 , N ，N, 分 别 是 它们 的 在 理想 条 件 下 的 最 大 容量 。 因 子 (1 - 一 ai | 反映 由 于 可 食 牧 草 对 环境 资源 的 
消耗 导致 的 对 它 本 身 增长 的 阻 河 作 用 和 毒 杂 草 对 环境 资源 的 消耗 对 可 食 牧 草 产生 的 影响 。 这 里 w 的 意义 
是 :单位 数量 的 毒 杂 草 ( 相对 WV, 而 言 ) 消 耗 的 供养 可 食 牧草 的 环境 资源 量 为 单位 数量 的 可 食 牧 革 ( 相对 WN, ) 消 
耗 的 供养 自身 的 环境 资源 量 的 w 信 。 a, 的 意义 是 :单位 数量 的 可 食 牧草 (相对 N, 而 言 ) 消耗 的 供养 毒 杂 草 的 
环境 资源 量 为 单位 数量 的 毒 杂 草 (相对 NN, ) 消耗 的 供养 自身 的 环境 资源 量 的 倍 。 模 型 (1) 的 4 个 平衡 点 


N(l-a Nll = a@; bs 衡 | 定 
1( 1) N,( 全 ] ,p,(0,0) ;这 4 个 平衡 点 的 稳定 性 分 析 结果 见 表 


要 
二 ot 1 = 而 问 5 


为 : pi(N ;0) Da(O5NV， ) » al 
1 中 证 明 过 程 见 附录 1) 。 


表 1 种 群 竞争 模型 的 平衡 点 及 稳定 性 


Table 1 Equilibrium point and stability of population competition model 


平衡 点 稳定 条 件 平衡 点 稳定 条 件 
Equilibrium point Stability condition Equilibrium point Stability condition 
pi (Ni,0) al < 1,a > 1 p> (0,N,) al > 1,a < 1 
N(l-a N,(1l-a 
| ne | a <1,a <l1 pa (0,0) 不 稳定 
1 -alas 1 — ala 
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2 毒 杂 草 入 侵 空间 扩散 模型 
毒 杂 草 和 可 食 牧 草 作为 两 个 物种 在 单个 斑 块 中 将 产生 竞争 ,而且 毒 杂 草 是 外 来 物种 ,在 生长 过 程 中 会 从 
密度 高 的 斑 块 向 密度 低 的 斑 块 发 生 迁 移 ,考虑 到 以 上 因素 建立 单个 斑 块 毒 杂 草 入 侵 空间 扩散 模型 
Wi _ Vi,j _ J 
di rixij( 1 六 | 六 1) 
(2) 
Ys ij Vij 
了 ri 人 1 TQ; N, 本 nN +A( Dy) 


模型 (2) 中 下 标 疙 表示 二 维 平面 上 的 网 格 中 的 第 ; 行 第 7 列 的 斑 块 ，( 也 yj, -4y;,) 为 第 忆 今 斑 块 中 
毒 杂 草 的 邻 体 对 第 疼 个 斑 块 的 入 侵 扩散 作用 ,本 模型 中 元 胞 自动 机 设计 采用 汉 诺 依 曼 型 :性 人 浴 元 胞 的 邻 体 只 
含 上 下 左右 4 个 邻 体 ,分 别 为 : W，( 左 邻 体 ) ,Ni ( 右 邻 体 )，N (上 邻 体 ) ,NA 下 邻 体 ) ,为 扩散 


系数 0 。 

2.1 具有 空间 扩散 的 毒 杂 章 入 侵 空间 扩散 模型 的 稳定 性 分 析 
讨论 模型 (2 ) 的 稳定 性 ,首先 要 计算 平 衔 点 , 令 ; 

Vi,j J 


-a 
N N 


x 外 


(xi = rixi(1 一 )=0 


Mi dn 
8 (Xij,)i)) 三 ri 人 1 一 02 一 “) 又 ( Sf 4y,,) =0 
N, N, 


则 平衡 点 为 : 
0 0, alN, alN, 
pi 0,5N, ;Pp2 0,5N, ; 也; 从. 三 A ;| N, — a o 
y 、 He /p= 0 Ry » 
定理 (证 明 过 程 见 用 录 2) 模型 (2 全 江 介 0, 人，| 稳定 的 条 作为 


0 >1,a, > ,ad <1; 


和 
0, 
平和 点 0, 和, 仆 确 巡 的 过 作为; 


2 
WE 


CN 
平衡 点 A 3 oa 稳定 的 条 件 为 : 


ai0 N,(a, - 1) 
ao > <1,aia, <min(2,— ”+1 ; 
N, 20, 
亚 逢 ; 占 aiN, 定 的 条 、 
平衡 点 py| N, - 040 稳定 的 条 件 为 : 
ai0s N, (a, 一 1) 
a >1,， <1,aa, <min{2,—————+1|。 
有 20, 
其 中 : 
4M 4M 
0=1+ /lt+ 全 二 = 二 
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下 
(a, — 1) + /0 -a,)’ -Ny (oo -1)M 


0, = 7 
NW (0 和 1) 


(oa -1) -= Ja -a,)’ - 一 一 (wa -1)M 
rN, 
0, = 5 
N (0 ey: 1) 


扩散 项 凡人 ya 屋 4y;,) 记 为 参数 M 。 

在 以 上 两 个 模型 中 ,模型 (1) 是 毒 杂 草 与 可 食 牧 草 的 种 间 竞 争 模 型 ,该 模型 只 考虑 了 当 两 个 种 群 在 同一 
个 自然 环境 中 生存 时 它们 的 相互 竞争 关系 。 而 模型 (2) 是 毒 杂 草 入 侵 的 空间 扩散 模型 ,该 模型 考虑 到 了 毒 杂 
草 和 可 食 牧 草 作 为 两 个 物种 在 单个 斑 块 中 将 产生 竞争 ,而 且 毒 杂 草 是 外 来 物种 仿生 生长 过 程 中 会 从 密度 高 的 
斑 块 向 密度 低 的 斑 块 发 生 迁 移 这 一 重要 因素 。 可 以 说 模型 (2) 是 模型 (1) 的 改进 。 时 得 到 的 结果 中 :没有 加 
和 人 扩散 项 时 ,模型 (1) 的 四 个 平衡 点 包括 零 平 衡 点 和 毒 休 草 灭 绝 的 平衡 点 以 及 一 个 共存 的 非 零 平 衡 点 ;加 入 
扩散 项 之 后 ,模型 的 零 平 衡 点 和 毒 杂 草 灭 绝 的 平衡 点 消失 ,而 且 毒 人 杂 草 与 可 丛 牧 章 的 共存 的 非 零 平 衔 点 由 一 
个 变 为 两 个 。 说 明 在 扩散 作用 下 ,两 种 群 的 可 能 出 现 两 种 稳定 情况 ,而 且 毒 杂 草 续 存 的 可 能 性 增加 。 


3 毒 杂 草 入 侵 的 空间 分 布 模拟 


将 单一 斑 块 中 的 毒 杂 草 与 可 食 牧 草 的 相互 竞争 扩展 到 空间 上 ,在 100 x 100 的 二 维 网 格 平面 模拟 两 物种 
( 毒 杂 草 与 可 食 牧 草 ) 的 种 群 动态 ,时 间 标 度 : 为 单位 时 间 步 长 。 初 始 时 刻 可 食 牧 草 与 毒 杂 草 的 数量 都 是 随机 
的 。 在 二 维 网 格 中 ,网 格 之 间 的 物种 会 发 生 迁 移 ,假定 物种 在 汉 诺 依 曼 邻 体 ( 即 只 考虑 相 邻 的 每 个 斑 块 边 的 4 
个 斑 块 ) 由 高 密度 向 低 密度 发 生 迁 移 ( 也 可 推 谍 汶 摩尔 邻 体 即 与 斑 块 边 和 顶点 相 邻 的 八 个 斑 块 发 生 迁 移 )。 
迁移 率 为 0.25 倍 的 密度 差 , 即 人 = 0.25。 

现 假定 m = 0.45 ,mr = 0.45 ,ai = 0.54o 30.7 ,最 大 容量 N 为 40, mV 为 60, 考 虑 到 毒 杂 草 为 外 来 物种 ,因此 
认为 该 物种 在 人 侵 过 程 中 具有 扩散 过程 而 作为 本 地 物种 的 可 食 牧 草 在 空间 上 是 没有 扩散 作用 的 ,通过 
Matlab 进行 编程 模拟 。 图 1 图 入 及 图 3 分 别 以 第 40,60,80 代 可 食 牧 草 与 毒 杂 草 空间 分 布 情况 为 例 给 出 空间 
第 40 并 毒 杂 草 空 间 分 布 第 40 代 可 食 牧 草 空间 分 布 


10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 100 


1 第 40 代 可 食 牧 草 与 毒 杂 草 空间 分 布 


Fig.l Spatial distribution of fortieth generations of edible grass and poisonous weeds 
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第 60 代 毒 杂 草 空间 分 布 第 60 代 可 食 牧草 空间 分 布 
100 pr pm 100 
90 上 90 
80 上 80 
70 70 
60 60 全 
o 
40 40 中 3 
30 30 mW 
20 上 光 20 
10 上 10 


80 90 100 


10 20 30 40 50 60 70 


2 第 60 代 可 食 牧 草 与 毒 杂 草 空间 分 布 


Fig.2 Spatial distribution of sixtieth generations of edible grass anid poisonous weeds 


第 80 代 毒 杂 草 空间 分 布 
wh 


3 第 80 代 可 食 牧草 与 毒 杂 草 空间 分 布 
Fig.3 Spatial distribution of eightieth generations of edible grass and poisonous weeds 


分 布 的 模拟 结果 ,其 中 白色 代表 生存 ,黑色 代表 已 灭绝 。 图 4 和 图 5 显示 :可 食 牧 划 生 物 量 大 概 在 18 代 之 后 
趋 守 不 变 , 而 毒 淋 草 的 生物 量 大 概 在 18 代 之 后 在 小 范围 内 波动 。 


4 毒 杂 草 入 侵 的 空间 分 布 假设 检验 


David & Moore(1954) 法 .Iwao(1968 ,1972) m-m 回归 法 ,平均 拥挤 度 指标 Lloyd，M.(1967) 等 方法 “六 是 
常见 的 种 群 的 检验 空间 分 布 状况 的 方法 ,根据 上 述 模拟 的 数据 ,利用 Matlab 编写 程序 ,判断 在 加 入 扩散 项 前 
后 毒 杂 草 的 空间 分 布 结 构 。 下 面 以 加 入 扩散 项 后 的 第 80 代 种 群 分 布 数据 为 例 进行 的 假设 检验 。 

(1) 利用 分 散 指标 C 对 数据 进行 检验 


基 : 沉 
根据 数据 计算 得 : 元 = 过 / SU 20.451834 


>» f 10000 
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可 食 牧 草 总 生物 量 


The total biomass of edible grass/ X 105 


0 1 1 1 


毒 杂 草 总 生物 量 
The total biomass of poisonous weeds/ X 104 


1 1 
0 10 20 30 40 50 
时 间 Time 


0 10 20 30 40 50 60 70 80 90 
时 间 Time 


4 可 食 牧 草 总 生物 量 随时 间 变 化 
Fig.4 The total biomass of edible grass changed with time 


又 因为 


所 以 
C 的 概率 为 95% 的 置信 区 间 为 


图 5 毒 杂 草 总 生物 量 随时 间 变 化 


Fig.S The total biomass offpoisonous, weeds changed with time 


>» (x; — x)? 


S = 148.05 
n-l 
S? (x -2¥ 
. 元 x.(n -3 
S? 
C = 一 =7123890084 
了 


2 
:2 | A NY 20 -1+0028287 
(n-W 99980001 


7.23890084 0 .028287 ,7.23890084 > 0.9717129 
显然 C 落 入 区 间 外 , 故 模拟 数据 属于 聚集 型 分 布 。 
(2) 利 用 David & Moore(195 物 礼法 对 数据 进行 检验 
首先 假设 模拟 数据 属 玉 随 机 分 布 , 则 有 


由 上 可 知 


2 2 2 
Moozs < 从 < Mo oo5 


[=—=C=7.23890084 


a | 部 


X=1(n -1) =7.23890084 x 9999 = 72381.7695 


所 以 


显然 不 接受 假设 ,所 以 模拟 数据 


i = 10000 - 1 = 9999 


Xo%s = 9723.732237 
Xaos = 10278.05632 


xX > Xo 0s 
属于 聚集 型 分 布 。 


2 和 
和 xX > X007s 
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(3) 利 用 平均 拥挤 度 指标 Lloyd，M. (1967 ) 对 数据 进行 检验 


n 


本 Zs > 。 545859.49 
平均 拥挤 度 指 标 参数 m= 二 一 一 一 一 ,判断 相应 数据 的 邻居 个 数 。 根 据 模拟 数据 求 得 m = a 
2 全 
26.69 , 即 平均 每 个 元 素 有 23.8569 个 邻居 。 因 为 m =x = 20.452 ,得 二 = 1.30500685 > 1, 所 以 模拟 数据 属于 


聚集 型 分 布 。 

(4) 利 用 Iwao(1968 ,1972) mm-m 回归 法 对 数据 进行 检验 

根据 Iwao 的 发 现 刀 二 Qw +Bm ， 因为 研究 的 是 个 体 分 布 , 所 以 分 布 基本 成 分 是 单个 个 体 ， 即 a=08 出 前 文 可 
知 m= 26.69, m =x=20.452, 所 以 可 得 6 ~ 1.30500685 > 1 ,很 容易 判断 模拟 数据 属于 聚集 型 分 布 。 

(5) 计算 Morisita's index 对 数据 进行 检验 


有 
由 公式 1 =n , 求 得 刀 = 1.30500685 ,自由 度 n -1 = 9999 ,HE 可 得 . 
(2*) - 之” 
Xo%s = 9723.732237 
X00s = 10278.05632 
, 计算 临界 点 
( Xo%ms — n+ D7 
C Uniform index ,三 = 0.99865 
- (ZZ) 7 
i Xo0ws Fn+ > sh 
ee Clumped index M. = = 1.0014 
和 (>, wx】 三 下 
所 以 
LL 三 M.>1 
它 计算 Standardized Morisita index 
6 L 到 M. 
= = 0.5 + 0.5 = 0.50001535 
(人 )》 站-= M. 


因为 五 = 0.50001535, 得 TI，> 0, 所 以 模拟 数据 属于 聚集 分 布 型 。 
采取 以 上 方法 ,通过 计算 机 模拟 数据 并 对 所 有 80 代数 据 进行 检验 ,得 到 的 结果 如 下 所 示 , 其 中 字母 “U” 
代表 均匀 分 布 , 字 母 人 ”代表 聚集 分 布 。 充 分 说 明 毒 杂 草 入 侵 可 食 牧草 时 ,以 较 强 的 苑 争 能 力 ,空间 分 布 型 从 
开始 的 均匀 和 从 布 逐渐 演变 成 最 终 的 聚集 分 布 。 
表 2 加 入 扩散 项 前 毒 杂 草 空间 分 布 结构 (1 到 80 代 均 为 A) 


Table 2 The spatial distribution structure of the poisonous weeds without diffusion (1th to 80th generations are A) 


代数 GenSrations 1 i 20 … 30 :… 40 1 50 … 60 … 70 … 80 
平均 拥挤 度 指数 Average congestion index A A Se A we A 
David & Moore 方法 David & Moore method A A A A A A A A 


A 


= 
= 


Morisita's 指数 法 Morisita's index 


由 表 2 可 以 看 出 ,在 加 入 扩散 项 之 前 , 即 扩散 系数 为 0 时, 毒 杂 草 在 第 1 代 到 第 80 代 一 直 为 聚集 分 布 。 

由 表 3 可 以 看 出 ,在 加 入 扩散 项 之 后 ,上 且 扩散 系数 取 0.25 时 , 毒 杂 草 在 第 1 代 至 第 24 代 时 为 均匀 分 布 ,从 
第 25 代 以 后 开始 变 为 聚集 分 布 。 

计算 1 至 80 代 加 入 扩散 项 的 植物 人 侵 模型 的 空间 分 布 相关 参数 , 毒 杂 草 从 第 1 代 开始 到 第 24 代为 均匀 
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分 布 ,从 第 25 代 开 始 其 空间 分 布 已 经 由 均匀 分 布 向 聚集 分 布 发 生 转变 ,此 后 一 直 是 聚集 分 布 。 综 合 两 个 表 可 
以 发 现 , 外 来 植物 物种 的 入 侵 及 空间 扩散 可 以 影响 植物 物种 的 空间 分 布 的 格局 ,减少 了 种 群 空间 分 布 的 聚集 
程度 。 


表 3 加 入 扩散 项 后 毒 杂 草 空间 分 布 结构 (1 到 24 代为 U,25 代 以 后 为 A) 


Table 3 The spatial distribution structure of the poisonous weeds with diffusion ( 1th to 24 generations are U, after the 25th generation are A) 


代数 Generations 1 :… 10 1 20 21 22 23 24 25 26 27 28 … 80 
平均 拥挤 度 指数 Average congestion index U U U U U U U A A A A I A 
David & Moore 方法 David & Moore method U 7 UU UU U U U A A A A :… A 
Morisita's 指数 法 Morisita's index U UN UUUUU A A A AIL… XK 


在 通过 Matlab 进行 编程 模拟 的 过 程 中 ,我 们 观察 到 ; 当 扩 散 系 数 、 固 有 增长 率 、 可 食 牧 革 沪 毒 杂 草 的 生存 
能 力 等 因素 变化 时 , 毒 杂 草 与 可 食 牧 革 的 空间 分 布 发 生 显著 的 变化 ,同时 , 毒 杂 草 与 可 食 牧 草 的 总 数量 随时 间 
的 变化 也 很 明显 。 当 取 扩 散 系 数 为 0, r = 0.45,r, = 0.45, a, = 0.5, a, = 0.7 时 ,由 于 没有 受到 扩散 因素 的 影 
响 , 毒 杂 草 从 第 1 代 开 始 到 第 80 代 都 是 聚集 分 布 ; 当 取 扩散 系数 为 0.1 ,其余 参数 不 变 时 ,由 于 受到 扩散 因素 
的 影响 ,可 食 牧草 与 毒 杂 草 的 分 布 结构 发 生 了 变化 , 毒 杂 草 在 第 1 代 至 第 14 八 时 为 狗 匀 分 布 , 从 第 15 代 以 后 
变 为 聚集 分 布 ; 当 取 扩散 系数 增加 到 0.2 时 , 毒 杂 草 从 第 一 代 至 第 /80 代 例 部 为 聚集 分布 ; 当 再 取 扩散 系数 为 
0.25 时 ,可 食 牧 草 在 第 1 至 24 代为 均匀 分 布 , 从 第 25 代 开 始 为 聚集 分 布 ;而 当 取 扩散 系数 继续 增加 到 0.3 时 ， 
由 于 扩散 系数 过 大 , 且 可 食 牧 草 与 毒 杂 草 竞 争 激烈 ,可 食 牧 草 在 很 短 的 时 间 内 灭绝 ,而 毒 杂 草 在 第 1 至 9 代为 
均匀 分 布 ,从 第 10 代 以 后 为 聚集 分 布 ; 当 扩 散 系 数 达 到 0.4 / 填 于 摘 散 系数 大 大 ,可 食 牧 草 很 快 就 灭绝 了 ,而 毒 
杂 草 在 第 1 代 时 为 均匀 分 布 ,从 第 2 代 开 始 便 属于 聚集 分 布 8 

从 以 上 的 模拟 过 程 中 我 们 观察 到 : 当 扩 散 系 数 从 0 逐渐 增 大 至 0.4 的 过 程 中 ,可 食 牧 革 的 灭绝 速度 越 来 
越 快 。 同 时 , 毒 杂 草 在 受到 扩散 因素 影响 的 时 候 , 其 均匀 分 布 保持 的 时 间 长 度 也 有 很 大 变化 , 当 扩 散 因素 从 0 
增加 到 0.25 时 , 毒 杂 草 均匀 分 布 保持 的 时 间 逐 渐变 长 ,而 当 扩 散 因素 从 0.25 增加 到 0.4 时 , 毒 杂 草 均匀 分 布 
保持 的 时 间 逐 渐变 短 。 


5 ”结果 讨论 


本 文 的 模拟 过 程 以 及 结果 充分 表明 ,在 生存 条 件 受 到 制约 时 ,种 群 的 发 展 与 其 生存 能 力 的 大 小 密切 相关 ， 
两 种 或 者 两 种 以 上 的 不 同 种 群 所 项 的 生存 资源 大 概 相同 时 ,在 同样 的 环境 下 ,生存 能 力 强 的 种 群 肯定 能 赢得 
种 群 竞争 的 胜利 5 。 就 如 同志 片 草场 中 的 牧草 和 毒 杂 草 , 它 们 的 生存 资源 基本 相同 , 毒 杂 草 的 生存 能 力 强 ， 
最 后 牧草 将 会 被 毒 杂 草 淘 汰 。 利 用 David & Moore( 1954) 方 法 .Iwao(1968,1972) m-m 回归 法 .平均 拥挤 度 指 
标 Lloyd ,M.(1967 ) 等 方法 对 毒 杂 草 的 空间 分 布 做 假设 检验 ,证 明了 毒 杂 草 的 人 侵 影响 了 草原 生态 分 布 , 毒 杂 
草 的 空间 分 布 型 也 南 最 初 的 均匀 分 布 演变 成 聚集 分 布 。 

本 文 在 假设 参数 时 ,认为 毒 杂 草 的 生存 能 力 略 大 于 牧草 , 故 在 模拟 中 参数 取 值 只 考虑 这 一 类 情况 。 由 于 
数据 是 模拟 出 来 的 ,可 能 并 不 能 很 贴切 的 描述 实际 的 情况 ,但 这 也 可 以 从 一 个 方面 来 了 解 研 究 毒 杂 草 空间 分 
布 ,进而 进行 有 关 的 毒 杂 草 问题 的 处 理 。 本 文 元 胞 自动 机 的 规则 是 自己 设 定 的 ,并 且 采 用 的 是 四 邻 域 , 这 个 是 
可 以 使 用 八 邻 域 .二 十 四 邻 域 , 得 到 进一步 推广 延伸 ,也 就 是 说 本 文 使 用 的 元 胞 自动 机 为 汉 详 依 曼 型 ,可 推广 
为 摩尔 型 以 及 扩展 摩尔 型 :9 。 
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ni =X 一 人 全 》 
附录 1 证 明 ;1) 对 pi(N,,0) 进行 平移 变换 |” ””“ ,代入 入 型 个 ) 得 : 


72 一 
dn riaiN, mA Ta 
qr = TN 一 N, ny Ww" N, Nn 
dn, ria,; 7 
dt = 疡 (1 JA) n, - Ny? 一 BW 
取得 线性 近似 系统 为 : 
dn riaiNi 
A Nn, 
dt N, 
dr 
站 交 =7,(1 -a,)n, 
riaiN, 
i 人 Wn Ty Ai =—7 
其 特征 方程 为 : N, = 0, 解 得 : 
A, =r,(1 -a,) 


0 A-r,(l-a,) 
当 ay < ja, 这 1 时 ,特征 根 有 负 实 部 , mm(N, ,0) 稳定 ,否则 不 稳定 (其 中 a4，< 1 是 根据 相 轨 线 分 析 的 结 
果 添 加 的 "这 也 与 生态 背景 相符 合 , 当 a，< 1 时 , 毒 杂 草 的 竞争 力 弱 于 可 食 牧 草 , 从 而 灭绝 ) 。 
9) 同 理 可 得 1, 当 a，> l,ao < 1 时 , p,(0,NN,) 稳定 ,否则 不 稳定 。 


Ni(l-a) 
NG-a) NI- WT 
3) 对 ps ， 进行 平移 变换 ,代入 模型 (1 ) 得 : 
1 一 ots 1 = gry (1 -a,) 
n> 一 一 1 二 

C102 

dz ma 一 1) riaNi(l -ai) mn 701 

去 Wi nN nn 
dt 1 -aa; N,(1 -aia,) Ni N, 
dn, ras Ns(1 -a,) r,(a, — 1) 72 _ 7202 


di Nilo) T1200 We NW? 


http://www.ecologica.cn 


